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We investigate the escape rate of a biaxial spin particle with an arbitrarily directed 
magnetic field in the easy plane, described by Hamiltonian Ti = —ASl — BS"^ — 
HxSx — HzSz, {A > B > 0). We derive an effective particle potential by using 
the method of particle mapping. With the help of the criterion for the presence 
of a first-order quantum-classical transition of the escape rate we obtained various 
phase boundary curves depending on the anisotropy parameter b = B/A and the 
field parameters ax,z = Hx^z/AS: azc(fec)'s, aa;c(fec)'s, and a^c = olzc{olxc)- It is found 
from a2c(&c)'s and axd^c^s that the first-order region decreases as h and ax (or a^) 
increase. The phase boundary line Uzc = otzc{otxc) shows that compared with the 
uniaxial system, both the first- and second-order regions are diminished due to the 
transverse anisotropy. Moreover, it is observed that, in the limit axe — > 0, a^c does 
not coincides with the coercive field line, which yields more reduction in the first-order 
region. We have also computed the crossover temperatures at the phase boundary : 
Tc{hc), Tc{axc,azc)- 

PACS number (s) :75.45.+j, 75.50.Tt 



1 



I. INTRODUCTION 



Recently there have been intensive studies on the quantum-classical phase transition of 
the escape rate in a single domain magnetic particle with many spins. In such a particle 
the magnetization direction of a collection of spins is oriented such that the magnetocrys- 
talline anisotropy energy is at a stable or metastable minimum depending on the existence 
of an external magnetic field. The escape from a stable or metastable state is governed by 
classical thermal activation rate, which is proportional to exp(— E'q/^b^)) |1| high tem- 
perature, and by quantum tunneling at temperatures below the energy barrier Eq. When 
these two escape rates are same there exists a crossover temperature Tq at which the transi- 
tion between classical and quantum regimes occurs. The transition can be either first-order 
or second-order. In the first-order transition the escape rate abruptly changes from the 
temperature-dependent thermal activation process to a practically temperature-independent 
quantum tunneling, so that the first derivative of the escape rate at the crossover temperature 
changes discontinuously. In the second-order transition, however, the escape rate changes 
smoothly from classical regime to temperature-dependent quantum tunneling (thermally as- 
sisted tunneling), and has a discontinuity of the second derivative. The determination of the 
transition order is closely related to the shape of the potential barrier which is controlled by 
the anisotropy constant and the external magnetic field. 

The quantum-classical transition of the escape rate was investigated by Affleck p[ and 
Larkin and Ovchinnikov. By using instanton technique they demonstrated that under 
certain assumption on the shape of potential barrier a smooth interpolation between the 
periods of oscillations at the bottom and the top of Euclidean potential well can be made, 
which leads to a second-order phase transition between classical and quantum regimes. 
Later, Chudnovsky 0] observed that the order of the phase transition in the crossover from 
classical activation to thermally assisted tunneling completely depends on the shape of the 
potential barrier. He has shown that the behavior of the energy- dependent period of oscilla- 
tions t{E) in Euclidean potential determines the order of the quantum-classical transition; if 
the Euclidean period increases monotonically with the increasing energy E from the bottom 
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of the Euclidean potential the transition is second-order, but if t{E) is nonmonotonic such 
that it has a minimum at E = Ei which is smaller than the potential barrier Eq the first- 
order transition occurs. More recently, Gorokhov and Blatter have obtained a sufficient 
condition for the first-order quantum-classical phase transition by looking at the behavior 
of the oscillation period in Euclidean time as a function of oscillation amplitude near the 
barrier top in the two-dimensional string model. In this case the first-order transition ap- 
pears when the amplitude- dependent period T(a), where a is the amplitude, is smaller than 
the zero amplitude period r(0) near the barrier top. This method has been subsequently 
extended to a quantum-mechanical model where mass has coordinate-dependence. ||^ 

The above approaches have been applied to the quantum-classical phase transition of 
the escape rate in a single domain spin system. Up to now two types of spin systems 
have been studied intensively: uniaxial and biaxial systems. For the uniaxial system, such 
as high-spin molecular magnet Mni2Ac, two models have been considered: one with a 
transverse field and the other with an arbitrarily directed field, described by Hamiltonians 
H = —DSl — HxSx and 7i = —DSl — H^Sx — H^Sz, respectively. The biaxial spin 
system, such as iron cluster Peg, [0] has attracted more attention, and several models have 
been taken into account. Liang et al. [|lT| considered a model without an applied field, 
Ti = K{Sl + ^Sy), (0 < A < 1) by using the periodic instanton approach and demonstrated 
that the coordinate-dependent effective mass plays an important role for the presence of the 
first-order transition. Based on the same approach Lee et al. |[T^ investigated the biaxial spin 



model with a transverse field, H = Ki^Sl + XSy) — HySy, and showed that the nonconstant 
mass which depends on both coordinate and field is important for the occurrence of the first- 
order transition. However, this approach involves a restricted range of applicability of A << 
1 for the biaxial spin system with a field. Such a restriction can be avoided by introducing 
the method of particle mapping. ||13[ The same model has been considered by Kim fl^ 



who used a quasiclassical method based on the particle mapping and found an analytical 
form of the phase boundary curve between first- and second-order transitions. The phase 
boundary of a biaxial system with longitudinal field described by 7-^ = K{S'^ + XSy) — H^S^ 
has also been obtained by Garanin and Chudnovsky |T5| employing a perturbation approach 



with respect to the transverse anisotropy, by Kim [|T4| using the quasiclassical method based 
on particle mapping, and by Park et al. [jl6| with the help of the Gorokhov and Blatter 's 
criterion and the particle mapping. 

In this paper we study the phase transition of the escape rate of a biaxial spin system 
with an arbitrarily directed magnetic field in the easy plane. In this case the transverse and 
longitudinal fields coexist, and there are three parameters which can be controlled by experi- 
ment: the anisotropy constant and two field parameters. In order to find a phase diagram for 
the transition orders these should be considered simultaneously. We will use the method of 
particle mapping to derive an effective particle potential from the spin Hamiltonian. Then, 
by applying the criterion developed in Refs. to this potential we will find various phase 
diagrams depending on the three parameters. Especially, we will obtain a phase boundary 
curve between the first- and the second-order transitions for the iron cluster Peg in which a 
relation between the transverse and longitudinal field parameters is given. We will also find 
the crossover temperature at the phase boundary. 

In the following section, we present a derivation of the effective particle potential based 
on the method of particle mapping. Here, we will take a constant mass so that the effective 
potential includes all the coordinate dependence. We will also review the sufficient condition 
for the presence of the first-order transition. In Sec. Ill we show various phase boundary 
curves between the first- and second-order transitions. Por completeness we will also briefly 
discuss two special cases: the biaxial systems with the transverse field only and with the 
longitudinal field only. These diagrams are compared with the previously obtained results in 
uniaxial system. The crossover temperature at the phase boundary is also computed here. 
Pinally, there will be a summary and discussions in Section IV. 



II. PARTICLE MAPPING AND THE CRITERION FOR THE FIRST-ORDER 

TRANSITION 

Consider a biaxial single domain magnetic particle with XOZ easy-plane and the easy 
Z-axis in the XZ plane. When an external field is applied along an arbitrary direction in 
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the XZ plane the Hamiltonian can be described by 



n 



'AS^ BS^ H-i-Srf HySy 



(2.1) 



where A and B are the longitudinal and transverse anisotropy constants, respectively satis- 
fying A > B > Our model is equivalent to 7i = K{Sl + \Sy) — H^Sx — HySy if we set 
A = K, B = [1 — X)K. For convenience, we introduce dimensionless transverse anisotropy 
parameter h = B/A{< 1) and field parameters = H^/SA, = Hz/SA where S is the 
spin number. Following Ref. this spin problem can be reduced to a particle moving in 
an effective potential. The equivalent particle Hamiltonian can be written as 

1 (P 



n 



2m dx"^ 



(2.2) 



where m = 1/2A is the particle mass, and V{x) = Av{x) is the effective particle potential 
given by 

. . 1 



vix] 



4dn X 



[S^ia^snx - a^cnx)^ - AbS{S + 1) 
— 2S'(2S' + l){bazSnx + a^jCnx)] 



(2.3) 



in which snx, cnx, dnx are the Jacobian elliptic functions with modulus k'^ = 1 — b. The 
Schrodinger-like equation corresponding to this Hamiltonian is 7^\l/(x) = i?\l/(x) where \E'(a;) 
is the particle wave function given by 



"^(x) = { - — I exp 



\dnx) 
X exp 



tanh ^(Vl — bsnx) 



tan 



2^6(1^ 



1 



1 — 6 cnx 
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(2.4) 



with 



$(x) 



s 

E 



/ snx + 1 



cnx 



(2.5) 



In Fig.l the effective particle potential is drawn. It has a metastable minimum with asym- 
metric barriers: small and large ones. The heights and widths of these barriers are governed 



by the anisotropy and field parameters 6, a^, ol^. For a given value of h tlie lieiglit of the 
small barrier decreases as or increases and vanishes at a critical value, while that of 
the large barrier increases. The critical value is determined by the metastability condition 



which, for the present model, is derived as ||T7| 



«l^ + «l^ = [2(l-&r^ (2.6) 

For a given value of h the metastable state exists inside the region closed by the curve 
(^xm = cxxnSs^zm)- We uotc here that because of the transverse anisotropy the metastability 
region of the biaxial system is smaller than the uniaxial system which corresponds to the 



case of 6 = 0. 11 



The escape rate at temperatures below the barrier height can be obtained by taking 
thermal average over tunneling probabilities. In the semiclassical approximation this can be 
expressed as 

/•-Bo 

Foe / dEP(E)e-^/^, (2.7) 
Jo 

where the metastable minimum is chosen to be zero energy. The tunneling probability P{E) 
can be approximated by the WKB exponent, P{E) ~ e""^*^^-* in which S{E) is the Euclidean 
action defined by 

rx2iE) I 

S{E) = 2V2^ / dxJv{x) - E, (2.8 

Jxj_{E) ^ 

where xi{E),X2{E) are the turning points corresponding to energy E. For large spin system 
the semiclassical approximation is well applicable so that we can neglect the contributions 
from paths fluctuating around the semiclassical saddle point path which minimizes the Eu- 
chdean action S{E). Within this approximation the escape rate becomes 

'S'Tmin(-£') 



exp ^ 

where STmin{E) is the minimum of an energy function 



(2.9) 



STiE) = SiE) + E/T, (2.10) 
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which is called thermon action or periodic instanton action in field theory. The condition 
for St{E) to have minimum requires dST{E) / dE = 0. Since the first derivative of an action 
with respect to energy in a potential well brings about the oscillation period we can write 

where negative sign attributes to the Euclidean potential in which the energy has negative 
values, and hence t{E) is an Euclidean time oscillation period. Since the Euclidean action 
S{E) is zero at E = Eq (i.e., at the top of the barrier) the minimum of the energy function 
becomes 

E 

STminiEo) = Sq = — (2-12) 

which is just the exponent of the thermal activation rate, i.e., the thermodynamic action. 0] 
The type of phase transition is determined by the behavior of t{E) with energy E. As has 
been well analyzed by Chudnovsky a first-order quantum-classical transition takes place 
when t{E) has a minimum at some energy Ei (< Eq). In this case the crossover temperature 
To = 1/t{Eq) is lower than Ti = 1/t{Ei), and there exists a temperature (Tq < < Ti) 
corresponding to an Euclidean oscillation period Tc at which the classical escape rate with 
5*0 and the quantum escape rate with the minimum thermon action STmm are connected, 
but their first derivatives with respect to T are discontinuous. That is, the semiclassical 
saddle point path jumps from the quantum regime to the classical one at temperature Tc. It 
can be noted from the behavior of t{E) that the first-order transition exist if the Euclidean 
oscillation period at an energy near the bottom of the Euclidean potential is smaller than 
r(-E'o), the period corresponding to the crossover temperature. 

The above argument can be generalized to the case in which the period is a function of the 
oscillation amplitude a near the bottom of the Euclidean potential. When E approaches to 
Eq the dynamics in the Euclidean potential well becomes small oscillations. At E = Eq the 
oscillation amplitude is zero, and the thermon action becomes the thermodynamic action. 
In this limit the period can be related to the crossover temperature To as following: 

To = l/r(0) = a;o/27r, (2.13) 



where r(0) is a period corresponding to zero amplitude, and Uq is defined as 



= f (2.14) 

V m 

where Xq is the position of the top of potential barrier. The solution of small oscillations near 
this point can be obtained from the Euclidean Euler-Lagrange equation, and the oscillation 
period of the solution can be expressed by r(a) = 27r/u;. From the above discussion the 
condition for the presence of the first-order quantum-classical transition is then given by 
r(a) < t(0). The difference r(a) — r(0) satisfying this condition can be calculated by 
the perturbation method in which the oscillation amplitude a is used as a perturbation 
parameter. |^ For the case of constant mass the condition becomes 

5 v"''^(xo) 1 ,,,,, , 

Below we will use this criterion to compute the phase boundary curve between the first- 
and the second-order transitions. Once we obtain the boundary curve we can also calculate 
the crossover temperature at the phase boundary from Eq.(13). 

III. PHASE BOUNDARY LINES BETWEEN THE FIRST- AND 
SECOND-ORDER TRANSITIONS 

We start this section by considering two special cases where either the longitudinal field 
or the transverse field is applied, from which we deduce some important features that can 
be directly applied to the case with arbitrarily directed field. We first consider the model 
with transverse field only. When the longitudinal field is zero the model can be described 

by 

n = -ASl-BSl-H,S,. (3.1) 

This is equivalent to the previously studied models; if we set A = K, B = {1 — \)K it is 
same as the model considered in Ref. ||12[, and ii A — B = K\\^ B = K± it becomes the 



model of Ref. WM. The present approach based on the criterion described above, however, 
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is different from theirs. In this case since = in Eq.(3) the effective particle potential is 
reduced as 

_ S^alsn'^x - 2S{2S + l)a^cna; - AbS{S + 1) 

Vt{x) — — , (3.2) 

4dn X 

and the corresponding wave function is given by 
^t(x) = ( 



cnx ^ ^ 



Vdnx 
X exp 



a^S -li b 1 

tan 



2^^! - b) 



1 — b cnx 



(3.3) 



We note that the effective particle mass has no coordinate dependence, which is different 
from previous models where the coordinate dependent mass played a crucial role to present 
the first-order transition. [0,113 In converting the spin problem to a particle one, whether 



the mass depends on the coordinate is a matter of how to set up the Schrodinger-like equation 
in the process of particle mapping (see Appendix). In our case all the coordinate dependence 
is included in the effective particle potential. As we can see below our approach gives the 
same results as Ref. 0]. 

The potential Vt{x) has now small and large barriers with same minima. Since the mass is 
independent of the coordinate it is obvious that the escape over the small barrier dominates. 
||T9| For large spin system, such as5'~5'+l~S' = 5' + l/2, the top of the small barrier is 
located at Xq = sn~^0 = 0. By equating both sides of the criterion Eq.(15) and evaluating 
the the derivatives of the potential Vt{x) at x = xq we obtain an equation of the phase 
boundary line between first- and second-order transitions: 



1 - 166c + 166^ + a/1 + 326c - 326^ 

atc(6c) = ; ^, , (3.4) 

^ ' 4(1 - 26c) 

where the subscript c represents that the anisotropy and field parameters are the values taken 
at the phase boundary. From Eq.(13), the transition temperature at the phase boundary Tc 
is then obtained to be 



'-tc 



1 



SA 27r\ 2(l-26c)' 



(3.5) 



where atdbc) is given in Eq.(19). These results are consistent with those in Ref. 0] if we 
reahze that = 2(1 — b)hx, b = kt/{l + kt). 

The phase diagram for the case with longitudinal field only based on the present approach 
has already been obtained by us, [|l^] so we just quote the results here: 



aiM = 2(1 - b. 



'I -2b, 
1 + b, 



Tic y/Sbc ll — be , . 

In Fig. 2 we have plotted the phase boundary curves atdpc), Oiidpc) and the coercive field 
line atm = otim = 2(1—6). The first-order transition exist below the curves ate (6c) and aic{bc)- 
In the region between these lines and the coercive field line the second-order transition lies. 
From these results we observe some important features. First, as the transverse anisotropy 
b increases both ate and aic decrease, and become zero at 6 = 1/2, i.e., the first-order region 
decreases in both cases. Thus, the biaxial system has smaller first-order region than the 
uniaxial system. Second, for a given value of b{< 1/2) the longitudinal field case has larger 
first-order region than the transverse field case. As we will see below these are common in 
the case with an arbitrarily directed field. 

The above features can be understood as follows. In the uniaxial system without trans- 
verse field the spin operator commutes with the Hamiltonian, and thus the spin becomes 
a constant of the motion (i.e., no dynamics). Moreover, in this case the barrier of the ef- 
fective particle potential becomes infinitely thick so that there is no tunneling, and hence 
no quantum-classical phase transition. In the presence of the transverse field, however, the 
uniaxial spin system becomes a dynamical one, and so the tunneling occurs. In the limit 
of very small transverse field the top of the effective potential barrier becomes flat, which 
is favorable to the first-order transition. When the transverse field increases, however, the 
situation becomes unfavorable to the presence of the first-order. Now, in the biaxial system 
the transverse anisotropy also gives dynamical origin to the spin system, and thus enforces 
the escape process from the metastable state, which leads to the suppression of the first-oder 
transition. Therefore, the decrease of azc(&c) is caused by the transverse anisotropy only, 
whereas the boundary line atc{bc) is affected by both the transverse field and transverse 
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anisotropy. 

We now consider the case with arbitrarily directed field. In this case the three parameters 
b, ax, oiz should be treated simultaneously, which is not a simple problem. In the present 
work we will fix one parameter and then compute the phase boundary with the other two 
parameters. We first calculate the phase boundary lines a2c(&c)'s for several values of a^-, 
and a2;c(&c)'s for different values of a^-, which are shown in Fig. 3. An immediate observation 
is that the first-order region for a given (or a^) diminishes as h increases, which shows 
the same trend as the = (or = 0) case. We also note that the first-order region 
becomes smaller for increasing values of a^^'s (or a^'s). This can be readily explained by 
the metastability condition in Eq.(6). When and b are given this condition yields the 
coercive field line for az such as azm{oix,b) = 2[(1 — 6)^/'^ — (a;^/2)^/^]^/^. From this it can 
be seen that the metastability region, where the escape process can be considered, is shrunk 
on the whole, which in turn suppresses the first-order transition region. 

Next, we compute the phase boundary for a fixed value of b. This is the case when a 
specific sample is prepared in experiment. Here, we take the molecular iron cluster Feg which 



has been studied in the previous experiment. [T^ In this case the transverse anisotropy 
parameter is given hj b = 0.29. In Fig. 4 we draw the phase boundary azdo^xc) and the 
metastability condition line C(zm{cixm) for b = 0.29. Since the metastability region of the 
biaxial system is smaller than the unixail system we can see both the first- and second-order 
regions become smaller than the uniaxial system. |^ The point of intersection with a^^-axis is 
larger than that of OzdO), which reveals that the case with longitudinal field only has larger 
first-order region than the transverse field only case . An interesting feature in this diagram 
is that unlike the uniaxial system azc does not coincides with the coercive field azm at zero 
transverse field. It can thus be realized that the first-order region is reduced more than 
the second-order region. In fact this is anticipated if we look into the Fig. 2 which displays 
that the longitudinal phase boundary line aic concurs with the coercive field line only at the 
point 6 = 0, i.e., the uniaxial system without transverse field. As is mentioned above this 
is the case of no tunneling, and thus we cannot think of any phase transition. Once b ^ 0, 
i.e., in the biaxial system, the two lines are separated such that the second-order transition 
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region becomes larger. 

The crossover temperature at the phase boundary is calculated numerically by using the 
values on the phase boundary lines in Fig.3 and the formula TcjSA = (l/27r) —v"{xq) /m. 
Fig. 5 shows Tzc/SA and Tr^c/SA as a function of be for several values of a^'s and a^'s, 
respectively. Since the be at a^c = decreases as increases (see Fig.3) the corresponding 
Tzc{bc)/SA ends at smaller value of be for increasing a^- The same trend can be found in 
Txcipc). It is noted that for a given be the crossover temperature at the phase boundary, 
Tzcipc) I SA^ rises as the transverse field grows, while T^Jyb^ / SA is lowered as the longitudinal 
field increases. By fixing the value of be and using Fig.4 we can investigate how the crossover 
temperature at the phase boundary varies with the field parameters. For be = 0.29 we have 
shown a 3d plot of Te{axe, o^zc) in Fig.6. From this picture it can be easily seen that Te is 
lowered with increasing aze, but rises as axe increases. 

IV. SUMMARY 

We have studied the phase transition of the escape rate of a biaxial spin particle with 
an arbitrarily directed magnetic field. By using the method of particle mapping we have 
obtained an effective particle potential to which the criterion for the presence of the first- 
order quantum-classical phase transition developed in Refs. has been applied. From 
this approach we have computed several phase boundary lines depending on the transverse 
anisotropy and field parameters. In the field vs transverse anisotropy plot, axe{be) or aze{bc), 
it is found that the first-order region decreases as be and ax^z increase. In the case of lon- 
gitudinal field vs transverse field plot, <yzc{<^xc), "we have observed that compared with the 
uniaxial system, both the first- and second-order transition regions are reduced (with more 
decrease of the first-order region) due to the transverse anisotropy. The crossover tempera- 
tures at the phase boundary corresponding to these diagrams have also been obtained. In 
the 3d plot of Te{axe, o^zc) we have observed that Te increases with axe, but decreases with 
azc- 

Experimentally, the phase diagram of azeic^xc) can be found in the octanuclear iron 
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cluster Feg in which S = 10 and 6c = 0.29. In this case the intersection with the a^j-axis 
in the phase diagram (see Fig.4) is given by a^c — 0.274 for which the magnetic field is 
estimated to be H^c = 0.64T, and for the longitudinal field, azdO) = 0.81 in the diagram 
for which Hzc = 1.9T, about three times of H^c- The corresponding temperatures are also 
estimated to be T^c — 0.39K which is the upper limit of the Tzc{bc — 0.29) 's for different 
values of ax''s, and T^c — 0.5K which determines the lower limit of T^dbc — 0.29) 's for 
different values of a^s. 



APPENDIX A: COORDINATE DEPENDENT MASS 

In converting the spin Hamiltonian of Eq.(l) to an effective particle Hamiltonian, if we 
set the Schrodinger-like equation to be 

1 (P 



+ V{x) 



(Al) 



2m{x) dx^ 

the coordinate dependent mass m{x), effective particle potential V{x), and the wave function 
^(x) are derived as following: 

1 



V{x) 



m{x) 



A 



2A(l + 6sinh2 x)' 



4(1 + 6sinh^x) 
X [alS^ - 4b{S^ + S-1/2)- 4a^S{S + 1/2) coshx 

+ 6^ sinh^ x-{b'^ + Ab{S^ + S-1)- alS^} sinh^ x 

— {4:bazS{S +1/2) coshx + 2axazS^} sinh x], 



*(x) = (l + ^sinh^x) 



X exp 
X exp 



CYzS _i / / r 

— tanh V 1 — otanh 



X 



tan 



1-b 



cosh a; 



$(a;) 



(A2) 



(A3) 



(A4) 
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with 

S (J 

= , = (sinhx + coshx)"^ ■ (A5) 

In this case, since the mass depends on coordinate the criterion for the first-order transition 
becomes more comphcated. Following Ref. the general first-order transition condition 
which includes the coordinate-dependent mass is given by 

[V"'{xo){gi + 92/2) + lV"%xo) + m'ixoHg^ 
+ m'{xo)u;^{gi + 92/2) + im"(xo)^o'] < 0, 



(A6) 



where 



_ oo^m'ixo) + V"'{xo) 

2m'ixo) + V"'ixo) 
- -A[Amixo)ui + V"ix,)y ^^'^ 

and Uq is the sphaleron oscillation defined as uJq = —V"{xo)/m{xo) from Eq.(15). Calculat- 
ing the derivatives of m(a;) and V{x) and substituting these into above condition we have 
same results as the constant mass case. 
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FIGURES 

FIG. 1. The effective particle potential with b = 0.3, = 0.2, = 0.3. The local minimum at 
X = —Xjn corresponds to a metastable state of the spin system described in Eq. (1). The inversion 
process of a spin magnetization vector takes place by the escape from the local minimum to global 
minimum at x = Xm along the path which passes through the small barrier. 

FIG. 2. The phase boundary lines atc{hc) and a/c(^c), and the coercive field line 
Oitm = aim = 2(1 - b). 

FIG. 3. The phase boundary lines (a) azcibcYs and (b) axcibcYs for several values of ax^s and 
a^'s, respectively. 

FIG. 4. Phase diagram for b = 0.29 which is the case of the iron cluster Feg. The solid line is 
the phase boundary azdctxc), and the dashed line corresponds to the line of metastability condition 
for b = 0.29. 

FIG. 5. Crossover temperatures at the phase boundary as a function of be- (a) Ti^/SA^s for 
ax = 0, 0.1, 0.3; (b) Ttc/SA's for = 0, 0.3, 0.7. 

FIG. 6. 3d plot of the crossover temperature Tc/SA at the phase boundary as a function of 
field parameters axe and a^c- 
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